Weak-Localization and Integrability in Ballistic Cavities 



Harold U. Baranger/b Rodolfo A. Jalabert/ 2 ^ 3 ) and A. Douglas Stone^ 
W AT&T Bell Laboratories 1D-230, 600 Mountain Ave., Murray Hill, NJ 07974 
(ty Physique de I'Etat Condense, CEA Saclay, 91191 Gif-sur-Yvette, France 
(3) Applied Physics, Yale University, New Haven, Connecticut 06520 

We demonstrate the existence of an interference contribution to the average mag- 
netoconductance, G(B), of ballistic cavities and use it to test the semiclassical theory 
of quantum billiards. G(B) is qualitatively different for chaotic and regular cavities, 
an effect explained semiclassically by the differing classical distribution of areas. 
The magnitude of G(B) is poorly explained by the semiclassical theory of coher- 
ent backscattering (elastic enhancement factor) — correlations beyond time-reversed 
pairs of trajectories must be included — but is in agreement with random matrix 
theory. 
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The main approach to relating the quantum properties and the classical mechanics of 
a system is semiclassical theory which expresses a quantum property in terms of the in- 
terference between certain classical paths [|TJ. The analytical semiclassical work to date 
includes only the interference between symmetry related classical paths — the diagonal 
approximation — while interference between paths unrelated by symmetry has remained 
largely intractable. In spite of this shortcoming, many results have been obtained; in the 
most studied example, the occurrence of chaotic or integrable classical dynamics determines 
the nature of the fluctuations in the quantum density of states [p]]. In the appropriate regime 
of many levels or short times, the result for chaotic systems agrees with random matrix the- 
ory, an approach to the statistical properties of complex quantum systems based on simple 
symmetry assumptions 

In this paper we demonstrate a new quantum interference effect — a change in the av- 
erage conductance of a ballistic cavity upon applying a magnetic field — connected to the 
breaking of time-reversal symmetry. Our semiclassical theory builds on those for studying 
quantum chaotic scattering [Q-Q and fluctuations in the quantum conductance 0. Within 
the diagonal approximation, we find that the average magnetoconductance is qualitatively 
different for chaotic and regular cavities. However, we find numerically that off-diagonal 
correlations — interference between classical paths unrelated by symmetry — make a large 
contribution to the magnitude of the average magnetoconductance and can in some cases 
eliminate it. In contrast, the random matrix theory of this effect does include off-diagonal 
correlations and therefore disagrees with the semiclassical diagonal-approximation result 
even though the time scale involved is rather short, the escape time from the cavity. Thus, 
this quantum interference effect highlights the successes and failures of the semiclassical 
approach currently used. 

The change in the average conductance upon applying a magnetic field is, of course, 
well known in disordered metallic conductors and is called weak- localization |7| . Despite 
strong similarities between ballistic chaotic systems and disordered systems, we show that 
weak-localization in ballistic systems is richer than in the diffusive regime: in addition to the 
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different behavior of chaotic and regular cavities, the effect of spatial symmetry and short 
non-ergodic paths is large. 

Since transport coefficients of microstructures directly measure scattering probabilities, 
they offer the possibility of direct experimental tests of "quantum chaos" |6|||. In fact, 
Marcus, et al. have recently reported || an experimental study of conductance fluctuations 
in microstructures in which a difference between nominally chaotic and regular shapes was 
observed. They also noted a large magnetoresistance peak at B = and connected this with 
weak-localization. 

We compute || the conductance of a cavity with two leads by relating it to the trans- 
mission intensity through G = (e 2 /h)T. Fig. 1 shows T(k) for a half-stadium structure. In 
order to reduce non-universal effects, we study an asymmetric structure in which a stopper 
blocks the directly transmitted paths; simpler structures are discussed below. The rapid 
fluctuations of T were studied previously [§,|n|; here we concentrate on the average con- 
ductance. A natural averaging procedure is to convolve with the derivative of the Fermi 
function to simulate non-zero temperature. The shift between the dashed and dotted curves 
in Fig. 1 is the average magnetoconductance that is the subject of this paper. 

The average T(B) in the inset of Fig. 1 shows that this ballistic weak-localization effect 
can be substantial. Traces for two structures are shown: the half-stadium structure and a 
similar structure with straight rather than curved sides. The classical dynamics in the half- 
stadium is chaotic while that in the straight-sided structure is regular. In the straight- 
sided structure the dynamics cannot be ergodic because the angle of a path exiting from the 
cavity must be related to the angle at entry through reflection from either vertical, horizontal, 
or diagonal walls. The difference between the behavior of these two structures — saturation 
versus linear increase — shows that ballistic weak-localization distinguishes between chaotic 
and regular classical dynamics. 

As a starting point for a semiclassical theory, we write T(k) in terms of classical paths 
which traverse the cavity PJT3|. For leads of width W which support iV modes, the total 
transmitted intensity summed over incoming (m) and outgoing (n) modes is 
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where 



F^W = V^A. exp [ifc(L fl - L u ) + i7r^,„]. (2) 

The only paths in the sum, labeled s and u, are those which enter at (x, y) with fixed 
angle sin 6* = ±m7r/kW and exit at (x',y ; ) with angle sin^' = ±rm/kW. In terms of the 
action S s , the phase factor is kL s = S s /h+ kysm.9— ky' sm9' plus an additional phase, 
4> s ,u, associated with singular points in the classical dynamics (see Ref. ||). The prefactor 
is A s = | (dy/d6') e ) \ /{WcosO'). 

Because the classical transmission coefficient is proportional to kW/n ||, we expect a 
linear contribution to the average quantum transmission and call the slope T. By averaging 



T \k) I (kW 1 7r) over all k |L4 |, one can show || that only terms with paired paths, s — u, 



contribute to T. The result, 
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f d(sin9) f d{sme') E A,, (3) 

8(0,9') 

is the classical probability of transmission. Thus, the leading order term in the average 
quantum conductance is the classical conductance ||. 

The quantum corrections are best discussed in terms of the reflection coefficient, R = 
N — T, for which identical semiclassical expressions hold in terms of reflected paths. The 
quantum corrections to R are 

= ^EE K% + E E (4) 

where we have separated the terms diagonal in mode number, Rd = J2n=i Rnn, from the off- 
diagonal terms. From results in disordered conductors, one expects that coherent backscat- 
tering will influence Rd{B) [/j. Previous work |l]-§J has shown that a typical diagonal 
reflection element is larger than a typical off-diagonal element by a factor of 2 when the 
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system is time-reversal invariant, a ratio known as the elastic enhancement factor. We pat- 
tern our discussion of 5R after previous semiclassical treatments of the elastic enhancement 
factor P,^|,P^| in which only the interference between symmetry related paths is included. 

Averaging 5Rd over all k ||14|| , (5Rd), we convert the sum over modes to an integral over 
angle, {it /kW)Y^ n ~^ I d6cos8, and note that the only fc-dependence is in the exponent 
so that all paths are eliminated except those for which L s = L u . In the absence of any 
symmetry, L s = L u only if s — u. However, time reversal symmetry at B = produces an 
additional contribution, namely by taking u to be s time-reversed. It is crucial in obtaining 
this pair to consider Ru in order that both paths satisfy the same boundary conditions on 
the angles. A weak magnetic field does not change the classical paths appreciably but does 
change the phase difference of the time-reversed paths by (S a — S u )/h = 2Q s B/(p where 
© s = 2it f s A ■ dl/B is the effective area enclosed by the path (times 2tt) and 0o — hc/e. 
Thus we arrive at the general expression 

(5R D (B)) = \f d( S m6) £ A s e l2e ° B ^ (5) 
1 s{e,e),s(e-e) 

which is a /c-independent contribution to the average quantum conductance. 

For a chaotic system, one can estimate both the magnitude and field scale of (5Rd)- 
First, if the mixing time for particles within the cavity is much shorter than the escape 
time, no preference is shown to scattering through any particular angle. To be precise, we 
assume the outgoing sin 6' are distributed uniformly and replace the sum over backscattered 
paths in Eq. ([5]) by an average over all sin#'. The resulting expression is the same as that 
for the slope of the average reflection, 71, in Eq. (|3|). Second, to estimate the field scale, 
we group the backscattered paths in Eq. (||) by their effective area and average over the 
distribution of this area, N(Q,8), 

/n/2 roc 
del dQN(Q,e)e i2@B/ ^. (6) 
-7r/2 J —oo 

Previous theoretical and numerical work has shown that in a chaotic system N(Q) oc 
exp (— « c /|@|) for large 6 independent of 9 where a c \ is the inverse of the typical area 
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enclosed by a classical path []16,^|,|6|,12|. Using this form for all G and the result for the 
magnitude, we find 

(6R D (B)) = K/[l + (2B/a d (f> ) 2 ]. (7) 

Note that the field scale can be much smaller than O through the area of the cavity. 

For a regular cavity, we estimate (5Rd{B)) by using the appropriate N(Q,6), following 
work in the energy-time domain For a fixed 9, we suppose that the trajectories 

are ergodic in real space and therefore that N(Q, 9) is exponential as in the chaotic case. 
However, unlike the chaotic case, the rate of decay depends on 9 and may vanish. For 
ergodic motion, this rate of decay is proportional to the square root of the typical escape 



rate 7 |16|j3|,p^] , N(Q,9) oc exp[— c|G| yj(9)], so that the points where 7^0 dominate the 
large behavior. Since in the regular structure of Fig. 1 7 vanishes linearly at 9 = ±7r/2 
as particles are injected close to a periodic orbit, one finds / d9N(Q,9) oc 1/0 2 . Taking 
the Fourier transform, we conclude that (5Ro{B)) oc \B\ for small \B\. (The unphysical 
cusp at B = is caused by deviations from 1/0 2 at very large 0.) Thus, a qualitative 
difference between chaotic and regular cavities results from the different classical distributions 
of effective areas. 

In Fig. 2 we compare the semiclassical predictions for the chaotic case to numerical 
results. The results for the structure with the stopper are in accord with the semiclassical 
theory: 5Rp is independent of k, its magnitude is within 25% of 71, and the elastic en- 
hancement factor is 1.97 ± 0.07 at B = and 0.99 ± 0.02 at B/a c i<po — 2. We emphasize 
that while the ^-independence of (SRd) is a general semiclassical prediction valid at large 
mode number, Eq. fl7|) depends on the separation of scales between the mixing time and 
the escape time, which is difficult to satisfy numerically. Some net variation as a function 
of k occurs in the structure with direct transmission paths (Fig. 2b) as well as a smaller 
total magnitude, indicating that short structure-specific paths can have a large effect on the 
average conductance. 

Fig. 2 shows that there is a large change in the off-diagonal reflection coefficients of 



opposite sign to (5Rd), a result not anticipated by the semiclassical theory above. Thus 
weak-localization is not equivalent to the coherent backs cattering (elastic enhancement) effect. 
This distinction does not appear to have been appreciated in much of the recent literature 
PJ^j9|,p!5[] . The importance of the off-diagonal contribution is even more apparent in the 
transmission coefficient since there are no time-reversal symmetric paths in the semiclassical 
expression for T. The off-diagonal term in Eq. (f|) is difficult to treat because the average 
on k does not eliminate interference of paths with a very small difference in action (extra 
factor of k from J2m) while the number of such paths is very large. More generally, analytic 
semiclassical calculations have only been able to treat interference between similar paths, 
and while this is adequate for aspects of the density of states |L] it is evidently inadequate 
for the magnitude of transport quantities. 

Despite this inaccuracy, the semiclassical theory presented here is successful in certain 
key respects. It shows that the reflection coefficients are sensitive to B through time-reversal 
symmetry, predicts the field scale correctly in terms of the average area enclosed by classical 
paths, and explains the difference between the chaotic and regular structures in terms of the 
distribution of the classical area. 

The semiclassical results suggest analyzing the numerical data by averaging the change 
in T(k) as shown in Fig. 3. The top panel demonstrates the difference between chaotic 
and regular structures: the curves for the half-stadia (chaotic) flatten out while that for 
the half-asymmetric-square (regular) increases linearly (except for very small B where it 
is quadratic). Not all of our chaotic structures show a clear saturation; however, all have 
rapidly changing magnetoconductance at small field followed by a more gradual rise. We 
attribute this deviation from the semiclassical theory to the small size of our structures. 

The results in the lower panel of Fig. 3 show a clear weak-localization effect for structures 
without stoppers. The error bars are larger than in the upper panel because of the greater 
variation with k produced by the direct paths. However, it is interesting, and important for 
experiments, that the direct paths do not mask the weak- localization effect: the difference 
between the chaotic and regular cavities is clear. The dotted curve in the lower panel shows 
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that the magnetoconductance is smaller for a symmetric stadium than for an asymmetric 
one, indicating the importance of spatial symmetries. 

Finally, the random matrix theory (RMT) for ballistic weak-localization [|TJ| proceeds 
by taking a chaotic quantum dot described by a Hamiltonian from the Gaussian Orthogonal 
Ensemble and coupling it to leads. For strong coupling, the result [jnj is (SR) = 0.25 while 
the elastic enhancement factor yields {SRd) = 0.5 = 1Z. The deviation of these results 
from those of the semiclassical diagonal-approximation (SDA) shows that this RMT does 
include off-diagonal contributions and is particularly interesting because the time scale, the 
escape time from the cavity, is smaller than naive limits on the validity of the SDA. If the 
coupling to the leads is weak, the RMT predicts a smaller weak-localization magnitude. We 
believe that the presence of short non-ergodic trajectories in our structures corresponds to a 
weakening of the coupling in the RMT. After including this effect, the predicted magnitude 
is 0.12 — 0.2 for our chaotic structures, in reasonable agreement with Fig. 3. 

In summary, the average magnetoconductance of ballistic systems shows great variety — 
the influence of chaotic versus regular classical dynamics, the effect of spatial symmetries, 
and the interplay between short and long paths. Understanding these effects presents a new 
challenge for the semiclassical theory. 

We thank B. J. van Wees for a stimulating conversation which helped initiate this work 
and appreciate discussions with P. Leboeuf, C. M. Marcus and A. M. Ozorio de Almeida. 
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8 



REFERENCES 

[1] For reviews see M. Gutzwiller, Chaos in Classical and Quantum Mechanics (Springer 
Verlag, New York, 1991) and M.-J. Giannoni, A. Voros, and J. Zinn- Justin, eds. Chaos 
and Quantum Physics (North- Holland, New York, 1991). 

[2] For a review of quantum chaotic scattering see U. Smilansky in Chaos and Quantum 
Physics, Ref. 0, pp. 371-441. 

[3] E. Doron, U. Smilansky, and A. Frenkel, Physica D 50, 367 (1991). 

[4] C. H. Lewenkopf and H. A. Weidenmiiller, Annals of Phys. 212, 53 (1991). 

[5] Y.-C. Lai, R. Bliimel, E. Ott, and C. Grebogi, Phys. Rev. Lett. 68, 3491 (1992). 

[6] R. A. Jalabert, H. U. Baranger and A. D. Stone, Phys. Rev. Lett. 65, 2442 (1990). 

[7] For reviews of weak-localization in metals see P. A. Lee and T. V. Ramakrishnan, Rev. 
Mod. Phys. 57, 287 (1985) and G. Bergmann, Phys. Rep. 107, 1 (1984). 

[8] C. M. Marcus, A. J. Rimberg, R. M. Westervelt, P. F. Hopkins, and A. C. Gossard, 
Phys. Rev. Lett. 69, 506 (1992). 

[9] H. U. Baranger, D. P. DiVincenzo, R. A. Jalabert, and A. D. Stone, Phys. Rev. B 44, 
10637 (1991). 

[10] R. B. S. Oakeshott and A. MacKinnon, Superlat. and Microstruc. 11, 145 (1992). 

[11] For reviews of classical chaotic scattering see T. Tel in Direction in Chaos, Vol. 3, ed. 
by Hao Bai Lin (World Scientific, Singapore, 1990) pp. 149-211 and Ref. @. 

[12] R. Jensen, Chaos 1, 101 (1991). 

[13] W. H. Miller, Adv. Chem. Phys. 25, 69 (1974). 

[14] It is important to average over an infinite window in k\ we use {A) = 
hm ? _ +00 (l / q) Jl c+q dkA(k) with q c W/n > 1. 

9 



[15] R. Bliimel and U. Smilansky, Phys. Rev. Lett. 69, 217 (1992). 

[16] M. V. Berry and M. Robnik, J. Phys. A 19, 649 (1986). 

[17] W. Bauer and G. F. Bertsch, Phys. Rev. Lett. 65, 2213 (1990). 

[18] S. Iida, H. A. Weidenmiiller, and J. A. Zuk, Phys. Rev. Lett. 64, 583 (1990) and Annals 
of Phys. 200, 219 (1990). 



10 



FIGURES 

FIG. 1. Transmission coefficient as a function of wavevector for the half-stadium structure 
shown in the bottom right. The T = fluctuations (solid) are eliminated by smoothing using a 
temperature (T = 6K for W = 0.5/xm) which corresponds to 20 correlation lengths. The offset of 
the resulting B = 2a c ;0o curve (dotted) from that for B = (dashed) demonstrates the average 
magnetoconductance effect. Inset: smoothed transmission coefficient as a function of the flux 
through the cavity (kW/n = 9.5) showing the difference between the chaotic (solid) and regular 
(dashed) structures. 

FIG. 2. Change in the total reflection coefficient (solid), as well as the diagonal (dashed) and 
off-diagonal (dotted) parts, upon changing B from to 2a c i<po. The curves are smoothed using 
a window of 1.5kW/n (30 correlation lengths for panel (a), 15 for (b)). The dashed ticks on the 
right mark the classical value of 1Z. Note the roughly ^-independent behavior of the curves in (a) 
and the large contribution of the off-diagonal reflection coefficients to the total weak-localization 
effect. 

FIG. 3. Weak-localization magnitude as a function of magnetic field for the six structures 
shown. The magnitude is obtained from (T(k, B) - T(k,B = 0)) k with kW/w £ [4, 11]. Note the 
difference between the chaotic and regular structures, as well as the sensitivity to symmetry in the 
lower panel. a c j is the inverse of the typical area enclosed by classical paths. 
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